The role of Joule heating in the formation of nanogaps by electromigration 
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We investigate the formation of nanogaps in gold wires due to electromigration. We show that 
the breaking process will not start until a local temperature of typically 400 K is reached by Joule 
heating. This value is rather independent of the temperature of the sample environment (4.2-295 
K). Furthermore, we demonstrate that the breaking dynamics can be controlled by minimizing the 
total series resistance of the system. In this way, the local temperature rise just before break down 
is limited and melting effects are prevented. Hence, electrodes with gaps < 2 nm are easily made, 
without the need of active feedback. For optimized samples, we observe quantized conductance 
steps prior the gap formation. 
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I. INTRODUCTION 

For single molecular electronics, electrodes are needed 
with a separation of typically one nanometer. In recent 
years, some creative methods have been developed to fab- 
ricate such devices. Some with tunable electrode dis- 
tance, like scanning probe methods ifilQand mechan- 
ically controllable break junctions, 0, 0, 13 and others 
with a fixed distance, like electromigration-induced tun- 
nel devices. 0, 0, HJ The latter have the advantage that 
they can easily be extended with a third (gate) electrode. 
Also, compared to other procedures, electromigration- 
induced nanogaps are easy to prepare, making this 
method one of the most promising techniques. 

Electromigration-induced nanogaps are produced by 
applying large current densities to gold wires. At high 
current densities j (typically 10 s A/ cm 2 ), momentum 
transfer from the electrons to the gold atoms leads to 
drift of the atoms, in the direction of the electron flow. 
This mass flux can lead to the growth of voids in the wire, 
finally leading to the formation of gaps. These gaps can 
have widths as small as a single nanometer, which makes 
them suitable for single molecular electronics. 0, Hi) 
However, to obtain gaps smaller than 2 nm, it is cru- 
cial that the process is indeed dominated by electromi- 
gration. Joule heating, resulting in melting and surface 
tension effects, can be the cause of much bigger gaps and 
gold island formation. 0, fl3L Il4| 

It is a challenging problem to prevent excessive heat- 
ing. The reason for this is that a high current density 
is only one of the requirements for nanogap formation. 
An additional condition is that the atomic mobility is 
high enough for substantial mass flux to occur. Since the 
mobility shows activated behavior, the local temperature 
plays a key role in electromigration. If the temperature of 



the sample environment is low (e.g., 77 K), Joule heating 
has to be substantial to induce the required local temper- 
ature. Below, we show, using a simple model and local 
temperature measurements, that gap formation typically 
takes place at a local temperature T start ~ 400 K. This 
value is rather independent of the temperature of the 
sample environment (4.2 K, 77 K, or 295 K). 
Unfortunately, the local temperature does not stay con- 
stant once the breaking process has started. While the 
slit is being formed, the current density and, hence, the 
local Joule heating increases. As a consequence, the 
local temperature diverges just before break down. In 
fact, temperatures can be reached up to the melting 
point of gold, leading to ill-defined junctions. This prob- 
lem can be solved by using an active feedback system, 
where the speed of the electromigration process is kept 
constant . [la. fl3t Il4 fl5L Il6fl In this way, the temperature 
is limited, giving better control on the final gap size. In 
this article, we discuss a simple model to describe the lo- 
cal temperature rise just before break down. We test this 
model by doing local temperature measurements during 
the electromigration process. We show that excessive 
heating can be prevented, by limiting the series resistance 
of the system. This makes a feedback system unneces- 
sary. By optimizing our sample design, we have increased 
the number of junctions with gaps < 2 nm from 15 to > 
90 percent. Furthermore, constrictions of only one or a 
few atoms can easily be achieved, leading to conductance 
steps of typically 2e 2 /h. 



II. THEORY: 



ELECTROMIGRATION IN 
METALS 
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Electromigration is commonly described as a mass flux 
under the influence of a high current density. However, 
a full description of the electromigration problem is not 
trivial. To give some insight, we refer to thermodynamics 
of irreversible processes. This theory considers all fluxes 
and forces involved, p^j In our case, there are three types 
of fluxes to be dealt with: the electron particle flux (J e ), 
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the flux of metal atoms (J m ), and the energy flux (J u )- 
These fluxes are induced by a set of three 'forces' (or 
potential gradients), Xy. For the particle forces, we can 
write Xj = — V/^ c . Here, p ec = jj, + Zetp is the elec- 
trochemical potential, with ip, p, and Z , the electrostatic 
potential, the chemical potential and particle charge (-1 
for electrons), respectively. The other 'force' is due to a 
temperature gradient: X u = V(l/T). From thermody- 
namics of irreversible processes we have following set of 
equations: 
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where the phenomenological constants Ly relate all 
fluxes to all forces, whilst obeying the Onsager relations, 
Lij — Lji. The coupled equations above provide a gen- 
eral description of the system. This description includes 
electromigration as well as thermo-electric effects and 
thermodiffusion. The latter refers to a mass flux due 
to temperature gradients (which can, in turn, be due to 
Joule heating) . Unfortunately, to solve the set of equa- 
tions for a three-dimensional geometry is a formidable 
task. With a few assumptions, however, a simplified rela- 
tion for the mass flux can be obtained. First, for materi- 
als with a high conductivity, we can ignore V/i e . Second, 
we choose to neglect thermodiffusion. In the experimen- 
tal part of this article, this assumption will be justified. 
Finally, we note that in practice all charge current is due 
to electron flux, so that — Vip — pj (with p the electrical 
resistivity). Hence, we obtain: 



J™ = -L*(Vp m - Z*e P3 ) 



(4) 



where we defined L* • = Li_j/T. Furthermore, we in- 
troduced an effective charge Z*: 



Z* = Z - 



J * 



T * 



(5) 



In words, the atoms behave as if they had a charge 
Z* . [13, [H, [l!j This effective charge is due to momentum 
transfer from electrons to atoms. Generally, >> Z, 
so the net force acting on the gold atoms will b"e in the 
direction of the electron flow. 



III. EXPERIMENTAL SETUP 

Two kinds of samples are fabricated on top of a 500 
nm thick SiC>2 layer on Si. For the first, called 'terrace 
samples', we use shadow evaporation and a resist bridge 
to obtain different thicknesses for the leads and the con- 
striction (see Fig. In the constriction, 15 nm Au 
is evaporated, while for the leads 150 nm Au on top of 
2 nm Cr is used. To decrease the constriction length, 
an extra 50 nm Au is evaporated. This way, a constric- 
tion of about 250 x lOOnm (length x width) is obtained. 
Due to thick gold leads, the total resistance is only 30 
Q (including the wires of the measurement setup). For 
the second kind of samples, called 'bow tie samples', we 
define a constriction in the shape of a bow tie with a min- 
imum width of 20 to 60 nm. In this case, the thickness 
of the gold in the leads and the constriction is the same, 
namely 17 nm. Furthermore, a 3 nm adhesion layer is 
used for most of the gold structure. By evaporating this 
layer at an angle, we make sure that the Chromium does 
not reach the constriction. Therefore, the very center of 
the constriction contains gold only. 

Electromigration is performed by applying a slowly in- 
creasing voltage (Keithley 230) to the wire, while moni- 
toring the current (Keithley 6517A). Measurements were 
done at room temperature (in air) and at both 77 and 
4.2 K (in cryogenic vacuum). After break down, tunnel 
currents were determined. In our setup, zero-bias tunnel 
resistances up to 10 12 f2 can be measured. 



IV. RESULTS AND DISCUSSION 
A. Excessive heating during electromigration 

A typical graph of an electromigration experiment is 
shown in Figure Here, breaking is performed for a 
'terrace sample' (Fig. ^l) by slowly increasing the voltage 
by 2 mV/s, at 77 K. The initial total resistance (including 
wiring) is 36 f2. At about 0.38 V, the current density 
has reached a value that causes substantial mass flux, 
starting the formation of a slit. At 0.4 V, the wire finally 
breaks down. Subsequently, a small tunnel current can 
be observed, indicating that the gap size is around 1 nm. 

A very different result is obtained for the experiment 
in Fig. |21 For this measurement, performed on a 'bow 
tie' sample, we added an extra series resistance. The to- 
tal initial resistance equals 1780 fi. Although the graph 
in Fig. Et) looks similar to the one in Fig^j), there 
are some obvious differences. First, the breaking pro- 
cess is much more abrupt (see insets). Second, the volt- 
age at break down is much higher (11 V vs. 0.4 V). 
The most dramatic difference is, however, that no tunnel 
current could be detected after break down. To inves- 
tigate the gap in more detail, we used atomic force mi- 
croscopy (AFM). In Figure |2h an AFM picture is shown 
of the 'bow tie' sample after the breaking process. The 
final gap size is approximately 100 nm, and the elec- 
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FIG. 1: a) Scanning electron micrograph of an electromigra- 
tion sample, 'terrace' type, made with electron beam lithog- 
raphy. The thickness of the gold layer is 15 nm (middle), 50 
nm, and 150 nm (leads). The constriction in the middle is 
about 100 nm in width, b) Representative breaking curve. 
The voltage is increased by 2 mV/s, while the current is ob- 
served. At 0.38 V, the electromigration process starts, leading 
to a decrease in current. At 0.4 V, after break down, a tunnel 
current can be measured of about 500 nA. Inset: magnifica- 
tion of region just before break down. Due to temperature 
increase and electromigration, dl/dV decreases slowly, until 
final break down occurs. Measurement at 77 Kelvin, the ini- 
tial total resistance is 36 Q. 
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FIG. 2: a) I(V) measurement during an electromigration ex- 
periment ('bow tie' sample). The total resistance before elec- 
tromigration is 1780 Q. At 11 V, the sample breaks down 
abruptly. No tunnel current could be observed afterwards. 
Inset: Magnification of region just before breaking. Measure- 
ment at 295 K. b) AFM picture of the sample after break- 
ing. As can be seen, the breaking mechanism for this sample 
was melting instead of electromigration. In between the elec- 
trodes, gold islands can be observed. 



trodes have clearly rounded off. Most likely this is due 
to local melting. Furthermore, there is a set of gold is- 
lands in between the electrodes. We note that the pres- 
ence of such nanoparticles may lead to Coulomb blockade 
effects. [lH HI l2fl EH] From Figure we conclude that 
the break down process has taken place in an uncon- 
trolled manner, leading to high local temperatures (and 
fields). We shall show below, that the breaking dynam- 
ics is strongly correlated with the resistance in series with 
the constriction, R s . We note that this resistance con- 
stitutes the main difference between the experiments in 
Figures Hand El Minimizing R s , results in a much better 
control of the final gap size, making a feedback system 
unnecessary. 

To demonstrate the significance of R s , we use a simple 
model for the constriction (a schematic drawing is shown 
in Fig. EJ). We consider a slit of length 1. Its width w(t) 
and height h(t) are decreasing functions of time, due to 
gap formation. The total resistance, Rtot, is the sum of 
the constriction resistance, R c (t) = j^^m , an d a series 
resistance R s . In a virtual experiment, we ramp up the 
voltage until electromigration begins at a bias V = V c 



(at t = 0). A more accurate way of saying that 'electro- 
migration begins' is to say that the flux of gold atoms, 
J m , reaches a certain critical value J^.[22] For this to 
happen, both the current density j and the mobility (re- 
lated to L* m m in eq. need to be substantial. The 
mobility is strongly affected by the local temperature, 
which increases due to Joule heating via p = pj 2 . Conse- 
quently, electromigration does not start, until a certain 
combination of current density j c and local temperature 
T s tart is reached. The latter is directly related to the lo- 
cal dissipation p c = pj 2 (below, we discuss this in more 
detail). 123 The voltage V c and local dissipation p c are 

related by p c = ^fk/(l + "r^oj) 2 - Next, we consider the 
breaking process for t > 0, keeping the voltage constant 
at V = V c - Due to electromigration a slit is formed in 
the constriction, leading to an increasing constriction re- 
sistance R c (t). Consequently, the local dissipation also 
increases, according to (in the diffusive regime): 
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FIG. 3: Schematic presentation of an electromigration sam- 
ple. The total resistance R consists of the constriction resis- 
tance R c and series resistance R 3 . Due to electromigration, 
the constriction will shrink in time, leading to a reduction of 
the width w and height h and an increase of R c . 



Let us first consider the limit with a small series re- 
sistance, << 1. In this case, the dissipated power 
is always equal to p c , independent of the width of the 
constriction. Hence, the local temperature stays close to 
T s tart and the gap grows slowly in time. In the other 
limit, j/^o) >> 1) the situation is different. At the start 
of the process (t=0), the local power equals p c . How- 
ever, the dissipated power increases rapidly during the 
breaking process, especially when R c (t) >> R c (0). In 
this limit, p >> p c . Consequently, uncontrollable heat- 
ing takes place. This can finally lead to melting of the 
electrodes. [2^ We note that equation JHJ) only holds in 
the diffusive regime, where the mean free path of the elec- 
trons is much smaller than the dimensions of the contact. 
Interestingly, as the electromigration process continues, 
the size of the constriction decreases. Provided that heat- 
ing is limited, one expects a cross-over from the diffusive 
to the ballistic transport regime. We will return to this 
phenomenon in section C. 

To have a better control of the final gap size, one has 
to limit the series resistance within the set-up as well 
as within the sample itself. For the latter we developed 
the 'terrace' samples that connect big leads to the con- 
striction (see Fig. [TJi). As compared to the 'bow tie' 
samples (see Fig. JSJd), this geometry has three advan- 
tages. The most important one is to limit the dissipated 
power during electromigration, by the reduction of R s . 
Another advantage is that the voltage needed to break 
the samples is limited. This is important, since, after 
break down all the voltage drop will be across the gap. 
This leads to a high electric field. It is well known that 
electric fields exceeding 2 V/nm, can cause a reorgani- 
zation of the electrodes. Depending on the shape of the 
electrodes, this reorganization can result in a larger gap 
size. [2jJ |2(| A third advantage of thick leads, is their use 
as heat sinks. In this way, the heat can easily flow away 
from the constriction. By the evaporation of 150 nm Au 
for the leads, we have a total resistance of typically 30 f2 
(including the wiring of the measurement setup) . 

To show the effect of a series resistance experimentally, 
a variable resistance is inserted in series with the sample. 
This way, we can check the final gap size as a function of 
R s . The result is shown in Fig. 0] We have investigated 
142 samples, of both types. The wires were broken by 
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FIG. 4: Percentage of samples with a zero-bias tunnel resis- 
tance < 10 12 Ohm after electromigration, as a function of the 
initial series resistance. High initial series resistances lead to 
excessive heating, causing large gaps. The number of sam- 
ples showing a tunnel current can be increased from 15 up to 
> 90 percent by limiting the series resistance. In total, 142 
samples were measured. The number of samples per bar is 
13/16/5/59/30/19. 

slowly increasing the voltage until a gap occurs. After- 
wards, the tunnel resistance at zero bias was measured, 
in order to make an estimation of the gap size. Zero bias 
resistances larger than 10 12 f2 could not be measured, and 
are counted as " no tunnelling" . As can be seen from Fig. 
2] there is a strong relation between the gap size and the 
total initial resistance. Considering samples with large 
series resistances, from 2 — 8 fcft, only 15 percent have 
nanogaps with a measurable tunnel current. Samples 
with low series resistances, however, give smaller gaps. 
Decreasing the value for the series resistance to about 30 
f2, gives nanogaps with a measurable tunnel current in 
more than 90 percent of the cases. From the tunnel cur- 
rent, an estimation can be made of the size of the gaps. 
These are in the range of to 2 nm. 



B. Local temperature during electromigration 

For the model discussed above, we assumed that gap 
formation does not start until a certain critical temper- 
ature T s tart is reached. To test if this is indeed the case, 
we performed measurements at 4.2, 77 and 295 Kelvin. 
We found no clear differences in final gap size. We did ob- 
serve, however, that a higher critical power p c is needed, 
when samples are broken at lower temperatures. This 
suggests that the local temperature at which electromi- 
gration is triggered, is more or less independent of the 
surrounding temperature. This is as expected, since the 
mobility of gold atoms is strongly dependent on the tem- 
perature (we will discuss this in more detail below). 

In order to find the local temperature during elec- 
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tromigration, T em , we fabricated four terminal devices. 
Knowledge of T em is important, since, for molecular elec- 
tronics, molecules are often put on the sample prior the 
breaking process. High temperatures could damage the 
molecules, right where the nanogap will be formed. 




5 10 15 
time [min] 



FIG. 5: Four-terminal resistance during electromigration. 
The voltage is set for an initial electromigration speed of 100 
(iQ/s, and kept constant afterwards. Due to electromigration, 
grain coalescence and heating, the resistance increases. Every 
couple of minutes, the resistance is measured at low voltage, 
and then V is set back to the initial value. From the change 
in resistance between low and high voltage, the average lo- 
cal temperature can be obtained (circles). The temperature 
of the environment is 77 K, while the temperature at which 
electromigration is initiated is 345 K. Just before break down, 
the temperature increases. The last temperature measured 
before break down is 450 K. The applied voltage is 0.72 V. 
The initial two-terminal resistance is 41 Q, the four-terminal 
resistance is 5.8 Q (at 77 K). The initial current density is 
ixW a A/cm 2 . 

To determine the local temperature, the samples were 
hrst cooled down to obtain the relation between tem- 
perature and (local) resistance (4.2 K < T < 295 
Kb[2?l l28l |29| The electromigration experiment is sub- 
sequently performed as follows. First, the voltage is in- 
creased till electromigration starts (dehned as dR/dt — 
100/ifi/s), and kept constant afterwards. An example 
is shown in Fig. [5] From here on, different processes 
will influence the local resistance. The most important 
ones are electromigration, heating and the coalescence of 
grains. To make an estimation of T em , it is crucial to 
deduce only the resistance change due to heating. This 
was done by interrupting the measurement every cou- 
ple of minutes, and determining the resistance at low 
voltage. |2^ Finally, the average local temperature T em 
is estimated from the resistance difference at low and high 
voltage and the (extrapolated) relation between temper- 
ature and resistance. 

For the trace in Fig. the initial four-terminal resis- 
tance was 5.8 (at 77 Kelvin). Setting the voltage to 
start the breaking process, increases the resistance dras- 
tically to about 8 (due to heating). Subsequently, the 
resistance increases in time, and accelerates just before 
breaking. Furthermore, the local temperature increases, 



TABLE I: Local temperature at which electromigration 
starts, T s tart, for 5 different samples. Interestingly, Tstart 
appears rather independent of the temperature of the sam- 
ple surroundings, To. For data set C and E, the temperature 
is also measured just before break-down, Tfinai- This tem- 
perature depends strongly on to what extent the process is 
controlled. Data set C is taken from Fig. |S]. 



Data set 


To[K] 


T s tart \K\ 


Tfinai [K] 


A 


4.2 


430 ± 40 




B 


4.2 


440 ± 40 




C 


77 


345 ± 10 


455 ± 10 


D 


295 


380 ± 40 




E 


295 


420 ± 10 


460 ± 10 



as predicted by eq. ©■ We note that there is always 
some series resistance, so that heating cannot be taken 
away completely. However, when high series resistances 
are used, the temperature does not only reach higher val- 
ues, it also diverges much faster. This is clearly shown 
in the insets of Figures and|2^- In Fig. QJj, a slow in- 
crease in differential resistance is observed. In contrast, 
for Fig. 03, gap formation goes all at once, and the 
breaking process cannot be resolved. 

In Fig. [5] the local temperature increases from 77 K 
to 345 K, when the voltage is set to start the break- 
ing process. The temperature determined just before 
break down, Tfi na i, is about 450 Kelvin. In Table jl], 
four other data sets are shown, at different environment 
temperatures. The temperatures at which electromigra- 
tion is triggered are close together, whereas the varia- 
tion in the temperature of the environment is about 290 
Kelvin. This is a strong indication that gap formation 
is indeed due to electromigration (at a critical tempera- 
ture) and not due to thermodiffusion via a temperature 
gradient. [29j 

To understand why T star t is rather independent of the 
temperature of the sample environment, To, we intro- 
duce a simple model. To fit its parameters, we use the 
data taken at 77 K (Fig. 0). We rewrite eq. ©, by 
noting that V/x m = Vtda/dx is the driving force due to 
stress cr: (19l |30I| 



- L *m,m{^v/ L - Z*epj) 



(7) 



Here Act is the stress built up over the length L. 
As long as Act < Aa ma x, the electromigration force is 
balanced by the stress gradient and the atom flux J m 
is zero. [13, till This defines a critical current density 
jmin — ^z'eph ' which is independent of temperature, 
since Z* p is temperature independent. 0, IH2 Using the 
Einstein relation D = L* m m dp/ dc, where D and c are 
the diffusion constant and concentration respectively, we 
can write down the temperature dependence of the atom 
flux. For the chemical potential we have /i = ferine. [T^ 
and for the diffusion constant D = D e~ Ea / kT so that we 
can write for the phenomenological coefficient: L* m m = 
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|yce _Ba / fcT . Here, E a is the activation energy of gold 
diffusion on the surface, which is 0.12 eV.|33 This leads 
to the following approximation for the atom flux due to 
the current density: 

Jm = ^(j-3min)e- E °' kT (8) 

where a = cD Z*ep/k < 0. For our measurements we 
use temperatures in between liquid Helium and Room 
temperature and the applied voltage over the constric- 
tion is typically 20 mV or higher. In this range, the 
inelastic scattering length of the electrons does not ex- 
ceed the contact size. Therefore, the effective tempera- 
ture in the contact can be described by T = T + [3j 2 , 
where j3 is a fit parameter.)^ We determine /3 from the 
first temperature point in Fig. [S] T start = 345 K at 
j = 4 x 10 s A/cm 2 . Furthermore, we estimate that j m in 
equals j min ~ 1 x 10 8 A /cm 2 . [l9| 

Equation [S] implies a strong relation between tem- 
perature and mass flux, especially due to the activated 
behavior. Gap formation starts as soon as the mass flux 
is large enough, i.e., when J m reaches a certain value J^. 
We define this quantity, in arbitrary units, as Jf n = 1. 
Going back to Fig. [SI with To = 77 K, we infer that for 
j = 4 x 10 8 A/cm 2 , we have J m = Jf n = 1. From this, 
we deduce the other parameter a. Knowing both a and 
P, we can plot J m as a function of current density j for 
all three experimental temperatures Tq, This is shown in 
Fig. For each To, we calculate the current density and 
temperature at which electromigration sets in, T stQrt , de- 
manding J m = 1. (For To = 77 K, this is by definition 
at j = 4 x 10 8 A/cm 2 with T = 345 K). For T = A.2K 
and T = 295K we find that T star t equals 330 and 420 
K, respectively. Hence, increasing the temperature of the 
environment by 290 K, increases T start by only 90 K. Re- 
lated to this is the fact that the critical current density 
increases with decreasing T . Although basic, the model 
is in rather good agreement with the values of Table [i] . 

Do previously applied molecules get damaged at these 
temperatures? As we have shown above, the local tem- 
perature just before break down depends on the ex- 
act sample design. Also, the maximum temperature 
molecules can stand, strongly depends on the molecule. 
At 460 Kelvin, which is Tfi na i in Table jl] , most molecules 
used in molecular electronics, do not get damaged. How- 
ever, care should clearly be taken. The electromigra- 
tion process should be controlled either by optimizing 
the sample design, or by using active feedback. If not, 
the local temperature increases dramatically as the wire 
gets thinner and may reach values up to the melting point 
of gold. 

C. Quantized conductance 

As we have shown, decreasing the series resistance 
gives a better control of the final gap size. Another in- 
teresting consequence of controlled electromigration, is 
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FIG. 6: Mass flux J m (arbitrary units) versus current den- 
sity j, for various surrounding temperatures To (see eq. JSJ). 
From Fig. 0, we know that electromigration starts (i.e., that 
Jm = 1) when j — 4 x 10 8 A/cm 2 with a local temperature 
Tstart— 345 K (rectangle). We use this to obtain the param- 
eters a and f3 in eq. (JSJl. With these, J m versus j, can also be 
plotted for To — 4.2 K and T = 295 K, respectively. Taking 
Joule heating into account (solid lines), the temperature at 
which electromigration starts is 330 K for To = 4.2 K and 
420 K for T = 295 K, respectively. Note that when Joule 
heating is neglected, electromigration can only set in at room 
temperature (dotted line). For 4.2 and 77 K, the critical mass 
flux (J m = 1) cannot be reached for reasonable j. 
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FIG. 7: Quantized conductance during the electromigration 
procedure at constant voltage. The voltage is set for an ini- 
tial electromigration speed of 100 fiQ/s, starting the breaking 
process. The resistance increases due to electromigration un- 
til, after 9 minutes, a conductance of 1 x 2e 2 /h is observed. 
After another minute, the configuration starts to reorganize, 
leading to a stepwise increase in conductance. Measurement 
at room temperature with V = 0.73 V and P= 1 x 10 -6 mbar. 



the possibility to observe the transition between the dif- 
fusive and the ballistic regime. If s °t formation is a 
well-defined process, at some point the constriction gets 
smaller than the inelastic scattering length. In this case, 
conduction becomes ballistic. This is indeed what we ob- 
serve in a number of cases. During the breaking process, 
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the sample conductance does not go smoothly to the tun- erlandse Organisatie voor Wetenschappelijk Onderzoek, 
neling regime, but locks in at plateaus equal to integer NWO, via a Pionier grant, 
values of the conductance quantum 2e 2 /h .(HE HEI An 
example is shown in Fig. \7\ Here, we start with a re- 
sistance of about 45 fl, increase the voltage until we see 
a change in conductance and, from there on, keep the 
voltage constant. Slowly, the conductance goes down, to 
finally stop at lx2e 2 /h. At this point, the electromi- 
gration process stops automatically, due to the almost 
perfect transmission of the gold atoms. The correspond- 
ing current density is about 57 /1^4/atom. This is close to 
the maximum value described elsewhere, |37] j and is much 
larger than the maximum current density in the diffu- 
sive regime. Nevertheless, the electrodes reorganize due 
to atomic diffusion. Slowly, the contact evolves from a 
single atom contact to a contact with 3, 5, 8 and 10 gold 
atoms. 38] At much higher voltages, it is still possible to 
break the constriction and obtain a tunnel resistance. We 
note that we have only observed quantized conductance 
in samples with minimized R s . If the power is unlimited, 
the high temperature prohibits the formation of a single 
atom contact. 



V. CONCLUSIONS 

To conclude, we have investigated the role of dissipa- 
tion during the formation of nanogaps due to electromi- 
gration. We find the following. First, some Joule heating 
is needed for gap formation to begin. The process does 
not start until a local temperature of typically 400 K is 
reached. This value is rather independent of the temper- 
ature of the surroundings. We relate this phenomenon 
to the activated behavior of the atomic mobility, which 
plays a crucial role in this diffusion process. Second, al- 
though Joule heating is important to start electromigra- 
tion, it can also lead to unwanted effects. If no measures 
are taken, the temperature in the constriction can in- 
crease up to the melting point of gold, during slit forma- 
tion. This leads to large gaps, possibly containing small 
gold islands. By minimizing the total series resistance 
of the system, we make sure the temperature keeps low 
during the electromigration process. In this way, we in- 
creased the number of samples with a gap < 2 nm to >90 
percent, without the use of active feedback. With the im- 
proved samples, gap formation is much slower. Hence, it 
is also possible to observe the transition from the diffu- 
sive to the ballistic regime. In a number of cases, contacts 
have been achieved that consist of only a single or a few 
atoms. 
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